Non-maximal 9 2 a, large #i 3 and tri-bimaximal 6 12 via 
quark-lepton complementarity at next-to-leading order 

Junpei Harada 

Niihama National College of Technology, Niihama 792-8580, Japan 
E-mail: haradaOsci . niihama- net .ac.jp 



ABSTRACT: We show that the next-to-leading order corrections in the quark-lepton comple- 
mentarity are important to explain the observed pattern of neutrino mixing. In particular, 
the next-to-leading order corrections 1) lead to a deviation of 023 from maximal mixing, 2) 
reduce the predicted value of sin 2 20i3 by 9.8%, 3) provide the same value of sin 2 6\2 as that 
of the tri-bimaximal mixing. This is shown by calculating sin 2 29ij(i, j = 1,2,3) to C(A 6 ) 
in the framework in which the product of the CKM and PMNS matrices is bimaximal. 
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1 Introduction and motivation 



The main recent developments on neutrino mixing[l-5] are related to the relatively large 

value of Bis, which was measured by recent experiments [6-10], and indications of significant 

deviation of 023 from maximal mixing[ll, 12]. In particular, non-maximal 023 is strongly 

Qh indicated by recent data, but the sign of 023 — tt/4 is not yet determined. Results of global 

(— | analysis are given in refs. [13-15]. 

To explain the observed pattern of neutrino mixing, the quark-lepton complementarity 
(QLC) has been widely investigated in the literature [16-41]. In particular, much attention 
has been paid to the class of models, £/pmns = ^ckm^' which can be obtained in grand 
CN unified theories (GUTs). The correlation matrix Vm is simply defined by the product of 

in 



the CKM and PMNS mixing matrices. In general Vm is not determined by theories. In 
this class of models, the observed two large mixing angles 0i2 and 023 indicate that Vm has 
two large mixing angles because all the mixing angles in Vckm are small. As the simplest 
possibility, for example, we can take Vm being the bimaximal mixing matrix Vb~ m , 

.^ UpmNS = VcKU V bm- ( L1 ) 

In this paper we only consider this minimal model. In ref. [25], expanding Vckm to 
C(A 4 ) (where A = sin0 c ~ 0.2253), we found that 

sin 2 20i3 = 2A 2 + 0(A 4 ) = 0.102 + 0(A 4 ), (1.2) 

sin 2 20 23 = l + O(A 4 ), (1.3) 

sin 2 20i2 = 1 - 2A 2 + 0(A 4 ) = 0.898 + 0(A 4 ). (1.4) 



These correspond to 0i3 ~ 0c/ V 2 — 9°, 023 — 45° and 0i2 — 36°, which are consistent with 
experimental results at leading order approximation. In particular, it is interesting that the 
predicted value of 1-3 mixing, 0i3 ~ 9°, has been confirmed by recent experiments [6-10]. 



However, the next-to-leading order corrections, the 0(A 4 ) terms in eqs.(1.2)-(1.4), 
should be calculated, because: 

• the 0(X ) corrections may not be enough small with respect to the recent experimental 
errors. Particularly, the Daya Bay Collaboration reported the precise value[9], 

sin 2 20i3 = 0.092 ± 0.016(stat) ± 0.005(syst) (la), (1.5) 

where the magnitude of systematic errors is the same order of that of A . 

• the same value of 1-3 mixing, #13 ~ 9°, can be obtained in various models which have 
different implications [12]. To distinguish models, we need a more precise prediction. 

• a deviation of #23 from maximal mixing is the 0(A 4 ) correction in eq.(1.3). To deter- 
mine the magnitude of deviation, the 0(A 4 ) corrections should be calculated. 

• the different values of 1-2 mixing have been predicted in various models[l-5, 11, 42]. 
In order to compare the various models in detail, we need more precise calculations. 

Motivated by these points, in this paper we calculate the C(A 4 ) corrections in the 
framework of eq.(l.l). We show that the next-to-leading order corrections 

• reduce the value of sin 2613 from 0.102 to 0.092, which is completely consistent with 
the result of Daya Bay Collaboration of eq.(1.5). 

• lead to a deviation of #23 from maximal mixing. The predicted value is sin 2 823 = 
0.446 + 0(A 6 ), which corresponds to #23 ~ 41.9°. 

• provide the value of sin 2 612 of 0.335, which is very close to the predicted value of 
tri-bimaximal (TBM) mixing of 0.333. 

The paper is organized as follows: in section 2 we present the explicit expressions 
of nine PMNS matrix elements to C(A 6 ). In section 3 we calculate the next-to-leading 
order corrections, and derive the neutrino mixing angles sin 2 29ij(i,j = 1,2,3) to C(A 6 ). 
Summary and conclusions are presented in section 4. 

2 PMNS matrix elements 

In this section we present the PMNS matrix elements U a i(a = e,p,,T,i = 1,2,3). We 
use the Wolfenstein parametrization[43, 44] of Vckm- The perturbative expansion by the 
parameter A = sin#c would make clear the consequences of eq. (1.1). 

The Vckm to 0(A 4 ) has been widely used in the literature, however, Vckm to 0(A 6 ) is 
necessary for calculations at next-to-leading order. After some calculations, we can write 
Vckm to 0(A 6 ) in terms of A, A, p and r/, 



/ l-f-¥ A A\* l + %)(p-ir]y 



Vckm 



2 8 v ' \j ' 2 

X + A 2 X 5 (l-p-irj) l-f-^(l + 4^ 2 ) AX 2 

V AX 3 (l-p-ir]) -AX 2 + AX\\-])-ir)) 1 - ^ 



2 

6 



+0{X b ), 

(2.1) 



where the Wolfenstein parameters A, A, p and r\ are defined by 

S12 = A, S03 = AX , sise = . . (2.2) 

y/T=J?[l-A*\4(j> + ifj)] V ; 

The quark mixing angles S{j = sin 0y and the Kobayashi-Maskawa CP phase 5 are defined 
by the standard parametrization. The CKM matrix written in terms of A, A, p and 77 is 
unitary to all orders in A. The values of the Wolfenstein parameters are given by[44] 

A = 0.22535 ±0.00065, A = 0.811±g:gf|, p = O.mjftoijj, fj = 0.345±g;gJ|. (2.3) 

The PMNS mixing matrix in the model of eq.(l.l) is 

(u el u e2 u e3 



t^PMNS 



U^U^U^] =vi KM V bm , (2.4) 

\u Tl u t2 u t3/ 



where the bimaximal mixing matrix V& m is defined by 




Vbm= I -3 g 4s I . (2.5) 



Substituting eq.(2.1) into (2.4), we can write the PMNS matrix elements to 0(A 6 ), 

1 A A 2 .4A 3 , , A 4 A 2 X 5 (\ \ mf ,^ . . 

A 4A 3 4 2 A 5 /I \ 

1,13 = ~vi + 7T (1 " ~ p + i!?) + ~VT (2 " p + 'y + c ' (A ' i) • (2 - 8) 

^=-i+^4G-^) + y{s + T +2i G-H} +o(A (2 - 9) 
^-s + ^-TG-^)-T{s + T +A G-H} +o(n (2l0) 



c, « = 71"7ll5 + - 4 J-72i8 + --' 4 l2- ? + i 'Jl + ( A6 »- < 2 '"> 

1 ^ 2 AA 3 . , .4 2 A 4 A\ 5 , . , ^,^6, . 010 , 

1 ^ 2 ,4A 3 , ^ 2 A 4 AA 5 , , „,, -, 

t/r2 = -2 + -y + "^-(P + ^) + "X" + ^7f (/f? + ^ +0(A } ' (2 ' 13) 

1 /1\2 42 \ 4 

The PMNS matrix given by eqs.(2.6) - (2.14) is unitary to all orders in A. 



3 Neutrino mixing angles 

3.1 013 

The recent experimental results on #13 from two accelerator experiments, T2K[6] and 
MINOS[7], and from three reactor experiments, Double-Chooz[8], Daya Bay[9] and RENO[10] 
were very important developments in neutrino physics. The global fit value of #13 is rather 
large, $13 ~ 9° for the normal and inverted mass hierarchy, and #13 = is now excluded at 
more than 10<r. 

The 1-3 mixing angle #13 is given by 

sin 2 2% = 4|C/ e3 | 2 (l-|^3| 2 ). (3.1) 

Eq.(3.1) indicates that #13 is determined only by U e3 . In eq.(3.1), #13 is a parameter defined 
by the standard parametrization, and the U e3 is the matrix element of any parametrization 
of t/pMNS [45] . Therefore, we can use the matrix elements given in section 2. 
From eq.(2.8), |C^ e 3| 2 is given by 

\U e3 \ 2 = ^-AX 4 (l-p) + 0(X 6 ). (3.2) 

Therefore, we find that 

sin 2 28 13 = 4|C/ e3 | 2 (1 - \U e3 \ 2 ) = 2A 2 - {1 + 4,4(1 - p)} A 4 + 0(X 6 ), 

= 0.1016- 0.0098 +C(A 6 ), 

A 2 term A 4 term 

= 0.092 + 0(A 6 ), (3.3) 

where we have used the values of A, A and p of eq.(2.3). Eq.(3.3) corresponds to #13 ~ 8.8° 
(an another solution of #13 ~ 81.2° is excluded, because sin 2 #13 = |C/e3| 2 < 0.5). 
The key points of eq.(3.3) are: 

• sin 20i3 is determined only by one parameter A at leading-order, and three parameters 
A, A and p are relevant at next-to-leading order. The r\ is irrelevant up to C?(A 6 ). 

• The next-to- leading order corrections reduce the value of sin 2 20i3 from 0.102 to 0.092, 
which is completely consistent with the result of Daya Bay collaboration of eq.(1.5). 
Thus, the value of sin 2Q\ 3 becomes smaller by 9.8% due to 0(\ ) correction terms. 

• The magnitude of the 0(X ) correction terms is approximately twice larger than 
that of the systematic errors of Daya Bay Collaboration, and therefore the 0(A 4 ) 
corrections cannot be negligible. 

• Because the C(A 6 ) terms are expected to be enough small, typically 0.1% corrections, 
eq.(3.3) would be tested by future experiments. Furthermore, eq.(3.3) can be used to 
distinguish various models which predict the same value of $13 at leading order. 
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Figure 1. Determination of the 1-3 mixing. Shown are the results from T2K[6], MINOS[7], 
Double Chooz[8], Daya Bay[9], RENO[10] and global fits of Forero et al.[13], Fogli et al.[14] and 
Gonzalez-Garcia et al.[15] for the normal hierarchy case (in the inverse hierarchy case the values do 
not differ by much). The QLC predictions are shown at leading and next-to- leading order. 

In figure 1, we summarize the theoretical and experimental values of sin 2 2(9i3. The 
results from T2K[6] and MINOS[7] are shown for sin 2 6 23 = 0.4 and 5 = it. For the QLC 
predictions, eq.(1.2) is shown for sin 2 2613 = 0.102 ± 5A 4 (leading order), and eq.(3.3) is 
shown for sin 2 2#i3 = 0.092 ± 5A 6 (next-to-leading order), which correspond to 0.089 < 
sin 2 2(9i3 < 0.114 and 0.091 < sin 2 26>i 3 < 0.093, respectively. Figure 1 shows that the QLC 
prediction at next-to-leading order is consistent with the results of five experiments [6-10] 
and three global fits[13— 15] at the la level. 

3.2 #23 

The recent global analysis indicates that there is a solid deviation of #23 from maximal 
mixing[ll, 12]. At present, the sign of 623 — vr/4 is not yet determined [13-1 5]. 

Eq.(1.3) shows that the 2-3 mixing is maximal at leading order. Eq.(1.3) also shows 
that a deviation from maximal mixing is not 0(A 2 ) but the 0(A) effect. Therefore, a 
deviation is not so large. To determine the magnitude of a deviation from maximal mixing, 
the 0(A 4 ) terms in eq.(1.3) should be calculated. 

The 2-3 mixing angle 623 is given by 



4|^ 3 | 2 |^ 12 
(l-|t4 3 P 2 ° 



sin 2 2c? 23 = " |L>3 ' U l 2 ■ (3-4) 



Eq.(3.4) shows that #23 is determined by three matrix elements, C/^3, U T s and U e ^. Of these 
elements, two are independent because of the unitarity condition |?7e3| 2 + |£^j3| 2 + |£A-3| 2 = 1. 
From eqs.(2.11), (2.14) and (3.2), \U^\ 2 , |C/ r3 | 2 and 1/(1 - \U e3 \ 2 ) 2 are given by 



IM 2 = \ ~ A 2 Q + A) + AA 4 (1 - p) + 0(A 6 ), (3.5) 

|[/ r3 | 2 = I +J 4A 2 + 0(A 6 ), (3.6) 

' 1 + A 2 + (^-2^1(1-p)1a 4 + 0(A 6 ). (3.7) 



(l-\Ue3\ 2 f ' I ^ 

Therefore, we find that 

sin 2 26> 23 = ^'^'2 = 1 - * Q + ^ 2 A 4 + 0(A 6 ), 
= 1- 0.0116 +C(A 6 ), 

A 4 term 

= 0.988 + C(A 6 ). (3.8) 

Eq.(3.8) indicates that the value of sin 2 2#23 is determined only by two parameters, A and 
A at next-to-leading order. The p and fj are irrelevant up to 0(A 6 ). 

To determine the sign of #23 — 7r /4, we calculate sin 2 #23 which is given by 

sin 2 23 = ., |C y, 9 . (3.9) 



l-|f/ e3 | 
From (3.2), 1/1 — |t/e 3 | 2 is given by 



' l + ^r + <!7-^(l-p)! > A 4 + 0(A ti ). (3.10) 



l-|C/ e3 | 2 2 [4 

Then, we find that 

** *■ = t^Cf = H* + ^ * - K* + *■) a * + o(A6) - 

= - - 0.0539 - 0.0005 +C(A 6 ), 

A 2 term A 4 term 

= 0.446 + C(A 6 ). (3.11) 

Eq.(3.11) indicates #23 < 7r/4. It is easy to confirm that eqs.(3.8) and (3.11) are consistent. 
The value of 6> 23 is #23 ^ 41.9°. 

In figure 2, we summarize the values of sin 023- The results of global analysis of Forero 
et al.[13] and Fogli et al.[14] are shown for the normal (NH) and inverted (IH) hierarchy 
cases. For the QLC predictions, eq.(1.3) is shown for 1 — 5A 4 < sin 2 2023 < 1 (leading order), 
and eq.(3.8) is shown for sin 2 2#2 3 = 0.988 ± 5A 6 (next-to-leading order), which correspond 
to 0.443 < sin 2 #23 < 0.500 and 0.444 < sin 2 #23 < 0.447, respectively. Figure 2 shows 
that the QLC prediction at next-to-leading order is consistent with the global fits[13-15] 
for #23 < 7r /4 at the 2<r level. 
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Figure 2. Determination of the 2-3 mixing. Shown are the results of global analysis from Forero 
et al.[13], Fogli et al.[14] and Gonzalez-Garcia et al.[15] for the normal and inverted hierarchy cases. 
The QLC predictions are shown at leading and next-to-leading order. 



3.3 6 12 

It has been known that a deviation of #12 from maximal mixing is large[l-5]. This can 
be naturally obtained in the model of eq.(l.l): eq.(1.4) indicates that a deviation of #12 
from maximal mixing is not 0(A 4 ) but the 0{\ 2 ) effect, and therefore a large deviation is 
obtained. For more precise prediction, in this subsection we calculate the C(A 4 ) terms in 
eq.(1.4). 

The 1-2 mixing angle is given by 



sin 2 20 



12 



A\TT |2|JT |2 
4|U e l| \Ue2\ 

a-\u e3 \ 2 ) 2 ' 

From eqs.(2.6) and (2.7), |C/ei| 2 and |C/ e 2| 2 are given by 

A 2 



U, 



el 



1 A 

77 + 



\u, 



V2 
A 5 

8^/2 
A A 2 



A ' f 1 AM s\ AX \l X 



(l + 2A){l + 2A(l-2p)} + 0(X 6 



'2 



V2 
A 5 



A ' f 1 am >\ AX \, ^ 



+ ^=(1 + 2A) {1 + 2A(l - 2p)} + 0(X b 



(3.12) 



(3.13) 



(3.14) 



Using these equations and eq.(3.7), we find that 

sin* 2^i2 = ^'~ e1 ' '^ ez ' = 1 - 2A 2 - 4A(1 - p)X 4 + 0(A 6 ), 



2nQ _ 4\U el ^\U e2 \ 2 
(l-|C/ e3 | 2 ^ 2 



= 1 - 0.1016 - 0.0073 +C(A 6 ), 

A 2 term A 4 term 

= 0.891 + C(A 6 ). (3.15) 

Eq.(3.15) shows that the value of sin 2 29\2 is determined by three parameters, A, A and p 
to 0(A 6 ). The fj is irrelevant up to 0(A 6 ). 

Eq.(3.15) corresponds to 812 — 35.4° (an another solution of #12 — 54.6° is excluded, 
because sin 612 = |C/ e 2| 2 /(l — | t/ e 3 1 2 ) < 0.5. The sin #12 written in terms of A, A and 
p is given in appendix A). It is interesting that the value of eq.(3.15) is very close to the 
tri-bimaximal (TBM) value of sin 2 20 12 = 0.889[46, 47]. 

In figure 3, we summarize the values of sin 2 #12. In addition to the results of global 
analysis[13-15], the predicted TBM value of l/3[46, 47], two proposed golden ratio (GR) 
values of 0.276(GR1[48, 49]) and 0.345(GR2[50, 51]) are shown. For the QLC predictions, 
eq.(1.4) is shown for sin 2 2#i2 = 0.898 ± 5A 4 (leaning order), and eq.(3.15) is shown for 
sin 2 2#i2 = 0.891±5A 6 (next-to-leading order), which correspond to 0.331 < sin 2 612 < 0.351 
and 0.334 < sin 2 612 < 0.336, respectively. Figure 3 shows that the QLC prediction at next- 
to-leading order is consistent with the global fits at the 2a level. 

4 Summary and conclusions 

In this paper we have calculated the next-to-leading order corrections in the framework in 
which the product between the CKM and PMNS mixing matrix is bimaximal, and have 
shown that these corrections were important to explain the observed pattern of neutrino 
mixing. Some conclusions obtained in this paper are: 

• 0i3 : the next-to-leading order corrections reduce the value of sin 2 28\% from 0.102 to 
0.092, which is consistent with five experiments, T2K[6], MINOS[7], Double Chooz[8], 
Daya Bay[9] and RENO[10], and global fits of Forero et al.[13], Fogli et al.[14], and 
Gonzalez-Garcia et al.[15] at the la level. A summary is shown in figure 1. 

• #23 : the next-to-leading corrections lead to a deviation of #23 from maximal mixing. 
A deviation of sin 2 2#23 from 1 is determined only by two Wolfenstein parameters A 
and A up to 0(A 6 ). The negative sign of #23 — 7r /4 is predicted, and the obtained 
value of sin 2 612 is 0.446, which is consistent with three global fits at the 2a level. A 
summary is shown in figure 2. 

• #12 : the predicted value of sin 2 812 at next-to-leading order is 0.335, which is very 
close to the TBM value of 0.333. These values are consistent with the global fits at 
the 2a level. A summary is shown in figure 3. 

Thus, the next-to-leading order corrections in the minimal model, eq.(l.l), are impor- 
tant to explain the observed pattern of neutrino mixing. It will be interesting to test the 
model by future precision neutrino experiments. 
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Figure 3. Determination of the 1-2 mixing. Shown are the results of global analysis [13- 15] and 
the predicted values of TBM[46, 47], GR1[48, 49] and GR2[50, 51]. The QLC predictions are shown 
at leading and next-to-leading order. 

A The 1-2 mixing and the matrix elements of the PMNS matrix 

For convenience, we present the sin #12, |t//xi| 2 , |^2| 2 , | C/-7-1 1 2 and |L r T 2| 2 in terms of A, A,p 
up to 0(A 6 ): 

sin 2 c? 12 = J- A_.A (1 _p )A 3 + ^Jl A{1 _p ) + A 2(^_p\\ +0{X 6 )i (A-1) 



2 a/2 V2 V2 18 v rJ \2 

,, 1 A A 2 (\ \ A 3 (\ \ AX 4 



^-I-T2 + ^\2 +A ) + T2\2- A )-— {1 - f) 

+ vH^ + T + <H} +0(n (A2) 

, r , ,2 1 A A 2 (\ \ A 3 (I \ AX* , , , 

-£{§ + T + A (H} + 0(A < A ' 3 > 

|2 _1 .4A 2 , .4A 3 p , AAWl ,A , _,.«, 



prtl =4- — + ^f + ^rU- A > +c,(A) ' (A ' 4) 



., 1 AA 2 A\ 3 p AAWl , 



p-i'-rT-^r-T r-'i^w (A ' 5) 
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